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Abstract 

The geometric phase associated with a many body ground state exhibits a signature of quantum 
phase transition. In this context, we have studied the behavior of the geometric phase during a 
hnear quench caused by a gradual turning off of the magnetic field interacting with a spin chain. 
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Geometric phases have been associated with a variety of condensed matter and sohd 
state phenomena [l-6| since its inception j^]. Besides, various theoretical investigations, 
geomteric phases have been experimentally tested in various cases, e.g. with photons js-lO]. 
with neutrons ll|, [l2| and with atoms 13|. The generation of a geometric phase (GP) is 
a witness of a singular point in the energy spectrum that arises in all non-trivial geometric 
evolutions. In this respect, the connection of geometric phase with quantum phase transition 
(QPT) has been explored very recently 14l4l6l|. There has been an increasing interest to 
investigate the importance of GP as a witness of QPT for various many body systems ItL 



231 ]. In essence, the small variations of external parameters resulting to the abrupt changes 



on the macroscopic behavior of a system is described by QPT. These critical changes are 
due to the presence of degeneracies in the energy spectrum and are characterized by long 
range quantum correlations. The geometric phase can be used as a tool to probe QPT in 
many body systems. 

Thermal phase transitions occur when the strength of the thermal fluctuations equals 
a certain threshold and the character of the stable phase changes. In contrast to that, 
zero temperature phase transitions, such as quantum phase transitions 2^, denote the 
crossover of different ground states at a certain critical value of some external parameter, 
where quantum fluctuations play a dominant role. Since response times typically diverge 
in the vicinity of the critical point, sweeping through the phase transition with a finite 
velocity leads to a breakdown of adiabatic condition and generate interesting dynamical 
(non-equilibriurn) effects. In the case of thermal phase transitions, the Kibble-Zurek (KZ) 



mechanism 
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, l26| explains the formation of defects via rapid cooling. This idea of defect 



formation in second order phase transition has been extended to zero temperature quantum 



phase transition (QPT) 27|, |28| in spin models. The KZ mechanism, explored in the one 



dimensional transverse Ising model helped to estimate the density of defects produced in 

_ 1 

QPT in the system|27l|. It is predicted that the density of kinks scales as ^, Tg being 
the quench time. For an adiabatic transition, i.e. for large Tg, smaller number of kinks are 
produced whereas for a fast transition, more number of defects are generated. 

From the existing literature, one may note that there is an increasing interest to in- 
vestigate the role of geometric phase in detecting QPT for many body systems from the 



geometrical perspective 



14j420l| and study the criticality of the system. The initial work 14 1 



showing that the geometric phase can be exploited as a tool to detect regions of criticality 
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without havingto undergo a quantum phase transition, was followed by subsequent works. 



It was shown ISjthat the geometric phase of the ground state obeys scaling behavior in the 
vicinity of a quantum phase transition and also the geometric phase can be considered as a 
topological test to reveal QPT [16]. In connection with QPT, there is a growing attention 



of studying the geoemtric phase of the ground state of various spin systems |17H20l|. Very 



recently, the study of geometric phase of the ground state for a comp 



icated inhomogeneous 
2l[ | there may exist more 



period-two anisotropic XY model in a transverse field showed that 
than one QPT point at some parameter regions and these transition points correspond to 
the divergence or extremum properties of the Berry curvature. 

These results motivated us to see how the geometric phase of the ground state of a spin 
chain behaves when the external magnetic field is varied. In this context, we have adopted 
the fascinating approach of studying the dynamics of a quench induced quantum phase 
transition, and obtain the dynamics of the geometric phase, which is marked as an indicator 
of QPT. Specifically, we have considered a spin chain with XY type of interaction in a slowly 
varying time dependent magnetic field and have studied the behavior of the geometric phase 
during a linear quench caused by a gradual turning off of the magnetic field interacting with 
a spin chain. In the adiabatic limit, we have derived the instantaneous geometric phases of 
the corresponding quantum states, and showed that the geometric phase is dependent on 
the quench time. 

We start with the review of the geometric phase of a spin 1/2 particle in a varying 
magnetic field. The Hamiltonian of a single spin 1/2 system in the presence of an external 
time dependent magnetic field is given by 

H{t) = ^B(t).a (1) 

where a = {ax,<Jy,crz) are the Pauli matrices and B(t) = Bon{9(t)) with the unit vector 
n = (sin 6' cos 0, sin ^ sin 0, cos 6*). If the external magnetic field is varied adiabatically the 
instantaneous energy eigenstates follow the directions of n and can be expressed as 



9 9 
t„;t> = cos-l tz> + e*'^sin-| 



9 9 
U]t> = sin-l t^> - e*'^cos-| 



(2) 



where | t2>) I 4-^> are the eigenstates of the operator. 



For a cyclic time evolution i.e. for B(T) = B(0), apart from the dynamical phase, the 
eigenstates acquire a geometric phase also and we can write, 



t„(T)>=e^V^^| tn(0)> (3) 



where the dynamical phase 6 = BQ{t)dt and the geometric phase (GP) "Jb = § ■ dX; 
A is the set of control parameters and A''^ = i <'|"„ | Va| tn> is the so called Berry connection. 
In our case, A = {6, 0) and a straight forward calculation shows that 

7t = -74 = 7r(l - cos6') (4) 

We should note here that though the eigenenergies depend on -Bo(^) the eigenstates depend 
on n(t) only. The GP is proportional to the solid angle subtended by B with respect to 
the degeneracy B = 0. It may also be noted that fluctuations in the external magnetic 
field will obviously induce fluctuation in the geometric phase (71- or 74,) of the corresponding 



spin through fluctuation in cos^ 29|]. The gradual slowing down of the magnetic field will 
induce an effect in the angle 6 which will subsequently affect the geometric phase of the 
corresponding state. In the present communication, we have adopted this idea in a spin 
system with XY type of interaction and investigated the dynamics of geometric phase when 
the system is evolved by a slowly varying time dependent magnetic field. 

Since XY model is exactly solvable and presents a rich structure, we extend our analysis 
for a chain of spin ^ particles with XY type of interactions to study the dynamics of geometric 
phase for an adiabatic evolution. This is a one dimensional model with nearest neighbor spin- 
spin interaction and the external magnetic field is allowed to orient along the z-direction. 
The Hamiltonian of this model is given by 

-AI ^ ^ 

where cxj's are the standard Pauli matrices at site i, a is the anisotropy parameter, = 
2M + 1 denotes the number of sites, B is the strength of the external magnetic field. We 
assume periodic boundary condition. 

uated by applying 



In this model, the geometric phase (GP) of the ground state is eva 
a rotation of around the Z-axis in a closed circuit to each spin [l^. 
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A new set of 



Hamiltonians if^ is constructed from the Hamiltonian ([5]) as 

= U{<P) H U\<P) (6) 



where 

+M 

U{4>) = n exp(i0a|/2) (7) 

j=-M 

and (7| is the 2; component of the standard Pauh matrix at site j. The family of Hamiltonians 
generated by varying has the same energy spectrum as the initial Hamiltonian and 
is TT-periodic in (f) With the help of standard Jordan- Wigner transformations, which makes 
the spins to one dimensional spinless fermions via the relation aj = {jX^^j (^t^ crj and then 
using the Fourier transforms of the fcrmionic operator 

dk — Ylij % 6xp ( ) with k — —M, ...+M the Hamiltonian H^j, can be diagonalized 
by transforming the fermionic operators in momentum space and then using Bogoliubov 
transformation. 

The ground state \g) of the system is expressed as 



1^) = Y[{cos^\0)k\0)-k - iexp(2i0) sin -j\l)k\l)-k (8) 

A;>0 

where |0)jk and |l)fc are the vacuum and single fermionic excitation of the k-th momentum 
mode respectively. The angle 9^ is given by 

cos k — B 

cos 9k = 7 (9) 

and 

Afe = yJ{cosk-By + a^ sin^A; (10) 

is the energy gap above the ground state. The ground state is a direct product of spins, 
each lying in the two-dimensional Hilbert space spanned by |0)jk|0)_fe and For 
each value of k, the state in each of the two dimensional Hilbert space can be represented 
as a Bloch vector with coordinates {2(f), 9k). The overall phase is given by the sum of the 
individual phases. The pseudomomenta k take half integer values: 

* = ±1H^, .±^^^ (11) 

The direct calculation shows that the geometric phase for the kth mode, which represents 
the area in the parameter space enclosed by the loop determined by (20, 9k) is given by 

rfc = 7r(l-cosefe) (12) 
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The geometric phase of the state \g) is given by 

r, = 7r(l-cos0fe) (13) 

k 

For an adiabatic evolution, if the initial state is an eigenstate, the evolved state remains 
in the eigenstate. So we may now derive the instantaneous geometric phases of this system 
due to a gradually decreasing magnetic field. Let us explore the situation when the system 
([5]) is driven adiabatically (slow transition ) by a time dependent magnetic field B{t) such 
that 

B{t < 0) = -— (14) 

B{t), driving the transition, is assumed to be linear with an adjustable time parameter Tq. 

Let the system be initially at time t(< 0) << Tg such that B{t) » 1. The instantaneous 
ground state at any instant t is given by 

\Mt)) = llicos^\0)k\0)-k-^exp{2^<|))sm^\l)k\l)-k (15) 

k 

We now use eqn. (12) and (9) to derive the geometric phase of the k*'^ mode which yields 

cos + — \ 
^k{t) = TT I 1 - (16) 
\cosk+^y + a^smH 



The variation of the geometric phase with time t, for a fixed k and a = 0.2 and for 
different values of the adjustable parameter Tq, the quench time, is shown in Figure 1. 

The geometric phase for an isotropic system with a = and quantum Ising model with 
a = 1 may now be easily obtained. 

For a = 0, rk(t) = 

(cos /c + — \ 
1 = (17) 

For a system of size the total geometric phase for the initial state is, 

^initial = Tfc (18) 

k 
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The magnetic field is gradually decreased by adjusting Tg and the critical point is attained 
at the instant t = —Tq with B = 1. Then the geometric phase for the k th mode is 

/ cos k + 1 \ , . 

Y ^/ [COS k + ly + sm"^ k J 
= for a = (20) 



TT ( 1 - -\/i(cosA;+ l)j for a = l (21) 

At the critical point the total geometric phase is 

^critical = ^^^k(t = —Tg) (22) 
k 

Finally, at t = 0, when the magnetic field is gradually turned off, the situation is a bit 
different. The configuration of the final state will depend on the number of kinks generated 
in the system due to phase transition at or near t = — and as such it will depend on the 
quench time Tg [26l] . The number of kinks is the number of quasiparticles excited at i? = 
and is given by 

AT = J^Pfc (23) 



k 



where Pk, t 
formula 31 1 



le excitation probability (for the slow transition) is given by the Landau Zener 



Pk ^ exp (-27rrgA;2) (24) 

As different pairs of quasiparticles {k, —k) evolve independently, for large values of r^, it 
is likely that only one pair of quasiparticles with momenta (fco, — ^o) will be excited where 
^o(= corresponds to the minimum value of the energy A^. Thus the condition for 
adiabatic transition in a finite chain is given by, 

» ^ (25) 
Hence, well in the adiabatic regime, the final state at t = is given by 

Q Q 

iV^/w) = |l)fco|0)-fco n (cosy|0)fc|0)_fc -2exp(2i0)siny|l)fc|l)_fc 

k,k^±kQ 

This state is similar to direct product of only — 1 spins oriented along {2(j),6k) where the 
state of the spin corresponding to momentum feg does not contribute to the geometric phase. 
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The total geometric phase of this state is given by 



rfinal{t^O)= Yl 7r(l-C0S^fe) (26) 



where 



cos k ,^„, 

^^"^'^ = / 2 ■ (27) 

V cos^ k + sm k 
Icos^fcl = 1 for q; = (28) 

|cos^fc| = COS A; for a — 1 (29) 

For a — 0, i.e. for an isotropic system, the geometric phase of the final state is given by 

^finalit ^0)^0 (30) 

and for a system with a — 1, i.e. for the well known quantum Ising model the total geometric 
phase of the state is given by 



Tfinalit^O)^ J2 7r(l-C0Sefc) (31) 

- 2 + cos ^) (32) 



2n ^ 



In a different case, if the system is driven slowly within the adiabatic limit and M number 
of defects are produced, the geometric phase of the final state will be given by 

r/i„„/(t = 0) = ^ 7r(l-cosA;) (33) 

k,k+kQ^,i=\,2,...M 

= 27r(A^ - 1) - 2A/'7r + tt ^cos ^ + sin ^ cot ^ - 1^ (34) 
If n be the average number of defects produced in the final state at 5 = 0, then 

N = nN (35) 

where n oc Hence, 

^finai{t = 0) = 27rN{n - 1) - 37r + tt (^cosTrn + sinTrA^cot ^) (36) 

Our analytical derivation shows that the in the adiabatic limit of a quench induced 
quantum phase transition, the geometric phase of the ground state in its final position with 
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vanishing magnetic field depends on the average number of kinks, and hence on the quench 
time. It is hoped that the present analysis may trigger investigations to explore the situation 
of physical realizations for obtaining the geometric phases associated with quench induced 
quantum phase transitions. 

To summarize, we have studied the dynamics of a quantum Ising model with a large 
number of spins and the geometric phase of the corresponding states. In our methodology we 
have adopted the novel idea of studying the dynamics of the geometric phase associated with 
a quantum phase transition when QPT is investigated in the light of KZ mechanism. The 
system initially at a high magnetic field was adiabatically evolved by gradually decreasing 
the magnetic field which is linearly adjustable by a time parameter r^, called the quench time. 
When the magnetic field is turned off at t = 0, the system, from its initial paramagnetic 
state, is finally in a ferromagnetic state with some defects generated in the system. The 
geometric phase of the final state depends on the average number of kinks produced in the 
system, which consequently depends on the quench time Tq. The present analysis traces 
the deformation of the geometric phase at criticality in a quench induced quantum phase 
transition associated with the generation of defects. 
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FIG. 1: Variation of with t for fixed a = 0.2 and k = tt/3 with different r^'s 
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